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Coastal  Subsea  Permafrost  and  Bedrock  Observations 
Using  dc  Resistivity 

PAUL  V.  SELLMANN,  ALLAN  J.  DELANEY  AND  STEVEN  A.  ARCONE 


INTRODUCTION 

Background 

Engineering  and  research  investigations  con¬ 
ducted  in  shallow  coastal  waters  often  depend  on 
geophysical  techniques  to  help  define  bottom 
conditions.  Seismic  techniques  have  commonly 
been  used  to  determine  layering,  changes  in  sedi¬ 
ment  type  and  depth  to  bedrock.  However,  there 
are  problems  with  noise  in  shallow  water  and 
there  is  high  attenuation  of  energy  where  sedi¬ 
ments  contain  gas.  Because  of  these  problems,  and 
the  lack  of  an  established  method  to  obtain  infor¬ 
mation  on  the  depth  to  the  top  of  shallow  ice- 
bonded  permafrost  in  Arctic  coastal  waters,  we 
initiated  this  study  using  dc  resistivity  techniques. 

Early  geophysical  investigations  of  permafrost 
used  dc  resistivity  methods  to  determine  its  distri¬ 
bution  on  land  (Joesting  1941,  Barnes  and  Mc¬ 
Carthy  1964).  More  recent  and  concurrent  devel¬ 
opment  work  (Scott  1975,  Corwin  1983,  Dyck  et  al. 
1983,  Sellmann  et  al.  1985,  Scott  and  Maxwell,  in 
press)  indicated  that  the  dc  resistivity  technique 
could  also  be  used  in  shallow  coastal  waters  (less 
than  30  m)  to  delineate  seabed  features.  The  dc 
method,  with  its  potential  for  good  resolution  at 
shallow  depths,  can  complement  deeper-sound¬ 
ing  Transient  Electromagnetic  Methods  (TEM) 
(Ehrenbard  et  al.  1983,  Walker  et  al.  1985).  Studies 
by  the  Earth  Technology  Corp.  (1985, 1986)  have 
suggested  that  modified  TEM  equipment  and 
procedures,  and  Magnetic  Induction  (MI),  might 
also  be  used  for  investigations  at  shallow  depths 
below  the  seabed.  At  present,  it  appears  that  a 
combination  of  dc  and  TEM  is  the  most  practical 
geophysical  approach  to  obtaining  information  on 
the  entire  offshore  permafrost  section. 

Researchers  interpret  resistivity  data  by  at¬ 
tempting  to  construct  layered  models  of  the  sub¬ 
surface  that  produce  resistivity  profiles  that  match 
field  observations.  Closed-form  inversion  solu¬ 
tions  are  virtually  impossible  to  achieve;  there¬ 
fore,  computational  programs  utilizing  fast  itera¬ 


tive  procedures  are  now  commonly  used  to  pro¬ 
duce  a  layered  model  to  fit  the  data.  These  pro¬ 
grams  require  estimates  for  the  resistivity  and 
thickness  of  each  assumed  layer.  Resistivity  data 
have  been  correlated  with  subsurface  sections 
having  a  known  distribution  of  ice-bonded  per¬ 
mafrost.  However,  much  of  this  information  is 
proprietary.  Near-shore  resistivity  data  are  also 
lacking  because  researchers  have  limited  experi¬ 
ence  with  shallow-penetrating  electrical  methods 
in  areas  of  ice-bonded  subsea  permafrost;  this  is 
one  of  the  reasons  for  this  study. 

The  electrical  properties  of  subsea  permafrost 
may  undergo  substantial  change  in  coastal  set¬ 
tings.  Inundation  can  dramatically  change  the 
properties  of  permafrost  by  exposing  it  to  salt 
water  and  warming  it  by  as  much  as  10°C.  Most 
subsea  permafrost  was  formed  on  the  continental 
shelves  when  they  were  exposed  during  periods 
of  low  sea  level  associated  with  times  of  major 
glacial  activity.  This  permafrost  was  then  covered 
by  rising  sea  level  as  the  glaciers  melted.  Even 
though  sea  level  is  now  fairly  stable,  permafrost  is 
still  being  inundated  in  areas  of  active  coastal 
erosion. 

Objectives  and  approach 

Our  objective  was  to  see  if  dc  techniques  could 
be  applied  to  seabed  studies  in  shallow  coastal 
waters  to  obtain  information  on  the  distribution  of 
bedrock  and  subsea  permafrost.  If  the  technique 
appeared  promising,  we  would  attempt  to  map 
the  top  of  ice-bonded  permafrost  and  establish 
resistivity  values  for  permafrost  in  areas  of  active 
coastal  erosion.  These  data  would  also  be  impor¬ 
tant  for  future  interpretation  of  electrically  based 
methods. 

The  initial  experimentation  was  done  at  coastal 
sites  in  New  England  to  obtain  field  data  in  sup¬ 
port  of  modeling,  thus  illustrating  the  influence  of 
water  depth  on  observations  made  with  a  floating 
cable  and  a  low  power  system.  The  results  encour¬ 
aged  us  to  construct  a  long,  floating  Wenner  array 


and  test  the  system  on  subsea  bedrock  features 
previously  mapped  by  seismic  methods  and  drill¬ 
ing  in  the  New  Haven,  Connecticut,  Harbor  (U.S. 
Army  Corps  of  Engineers  1981).  This  site  was  also 
selected  since  the  geologic  setting  (mud  over  bed¬ 
rock)  provided  a  good  electrical  approximation  of 
bottom  conditions  in  an  area  of  subsea  permafrost. 

Positive  results  at  the  New  Haven  site  led  to  the 
studies  near  Prudhoe  Bay,  where  ice-bonded  per¬ 
mafrost  is  known  to  exist  at  shallow  depths  below 
the  seabed.  Several  sites  with  different  coastal 
relief  and  beach  configuration  were  studied.  We 
assumed  that  the  general  distribution  of  subsur¬ 
face  materials  would  be  similar  at  all  study  sites 
and  that  coastal  erosion  rates  were  approximately 
1  m/yr.  We  collected  data  on  land  with  dc  sound¬ 
ings  made  parallel  to  the  shoreline,  using  elec¬ 
trodes  driven  into  the  ground.  Observations  ex¬ 
tended  progressively  inland  beyond  any  evidence 
of  flooding  of  the  tundra  surface.  We  varied  the 
spacings  between  soundings  on  land  to  provide 
coverage  of  local  relief  and  features  found  on  the 
beaches,  bluffs  and  adjacent  tundra  surfaces.  The 
farthest  inland  soundings  were  considered  refer¬ 
ences,  representing  subsurface  electrical  proper¬ 
ties  not  yet  modified  by  coastal  processes,  and 
were  made  on  what  appeared  to  be  undisturbed 
tundra  surfaces.  A  floating  cable  was  used  to  ex¬ 
tend  the  survey  lines  up  to  2.8  km  offshore  by 
profiling  and  sounding  at  regular  intervals  from 
the  beach.  The  cable  was  normal  to  the  shoreline 
for  most  of  this  surveying. 


DC  SOUNDING  METHOD 

Resistivity  was  measured  using  the  Wenner 
array  (Fig.  1 ),  which  consists  of  four  colinear,  equis- 
paced  electrodes.  The  two  outer  electrodes  inject 
and  receive  current  I,  and  the  resulting  gradient  in 
potential  V  is  measured  between  the  inner  elec- 


Figure  1.  Multielectrode  Wenner  array.  The  electrode 
separation  a  was  varied  from  0.5  to  50  m  on  land  and  3 
to  50  m  offshore. 


trodes.  For  homogeneous  ground,  the  real 
resistivity  p  in  ohm-meters  is  equal  to  the  meas¬ 
ured  apparent  resistivity  pa,  which  is  found  from 
the  formula 

pa  =  2rc  a  ¥ 

I 

wher ea  is  the  interelectrode  spacing.  Over  layered 
earth,  pa  is  also  a  function  of  the  layer  thicknesses 
and  resistivities;  therefore,  it  is  an  apparent  quan¬ 
tity.  Computational  modeling  must  then  be  used 
to  interpret  the  sounding  curve. 

Two  of  the  CRREL  cables  were  fabricated  from 
standard  single-conductor  wire.  Cable  1  was  210 
m  long  with  seven  a  spacings  that  ranged  from  10 
to  60  m.  It  was  used  for  both  bottom  and  floating 
cable  measurements.  Cable  2  was  a  short  bottom- 
cable,  30  m  long,  with  six  a  spacings  that  ranged 
from  2  to  11  m.  The  third  cable,  fabricated  by  a 
commercial  cabling  company,  had  a  polyurethane 
outer  jacket  and  was  bouyant.  It  was  150  m  long 
with  eight  Wenner  a  spacings  from  3  to  46  m.  The 
progressively  larger  a  spacings  in  the  Wenner 
array  telescoped  out  from  the  end  nearest  the  boat 
so  that  the  first  electrode  was  always  in  use.  All  of 
the  cables  used  simple  lead  electrodes  as  described 
by  von  Arx  (1962).  These  cables  provided  a  wide 
range  of  electrode  configurations  and  separations 
for  both  floating  and  seafloor  measurements. 

Modem  equipment  and  operating  techniques 
are  designed  to  overcome  two  major  obstacles  in 
making  dc  resistivity  measurements:  1)  excessive 
resistance  caused  by  poor  contact  between  the 
electrodes  and  the  ground,  which  can  occur  in  dry, 
coarse-grained  soil,  in  frozen  ground  or  in  rock, 
and  2)  measurement  interference  associated  with 
natural  earth  currents  and  potentials.  Contact 
resistance  is  not  a  problem  in  the  marine  environ¬ 
ment,  since  the  electrodes  are  immersed  in  seawa¬ 
ter.  On  land  we  used  copper-clad  electrodes  driven 
approximately  15  cm  into  the  ground  at  a  spacings 
up  to  50  m.  The  low  resistivity  of  the  thawed  active 
layer  and  marine  sediments  on  the  beach  ensured 
good  electrode  contact. 

A  regulated  current  signal  with  a  total  period  of 
32  seconds  was  generated  by  a  Huntec  transmitter, 
and  potential  was  measured  with  a  Fluke  digital 
voltmeter.  The  Huntec  unit  has  switch-selectable 
output  impedance,  and  output  current  is  variable 
to  1.5  A.  In  most  cases  the  maximum  output  cur¬ 
rent  was  used  for  the  marine  surveys.  An  illustra¬ 
tion  of  this  instrumentation  and  one  of  our  marine 
cables  is  shown  in  Figure  2. 
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Figure  2.  Generalized  illustration  of  multiconductor  cable,  electrode  geometry  and 
switching  equipment  used  for  the  offshore  observations.  Switching  equipment  allows 
selection  of  a  separation. 


COMPUTER  MODELING 
AND  INTERPRETATION 

Evaluation  and  interpretation  of  data  were 
based  on  the  well-developed  theory  for  resistivity 
sounding  of  layered  media.  We  did  the  forward 
modeling  used  to  predict  the  results  of  the  earlier 
studies  at  the  New  England  sites  with  the  aid  of 
theory  developed  by  Mooney  et  al.  (1966).  The 
Beaufort  Sea  data  for  the  later  studies  were  in¬ 
verted  using  a  commercially  available  (but  propri¬ 
etary)  program  called  RESIX  (Interpex  Limited, 
Golden,  Colo.).  Such  programs  have  evolved  from 
the  work  of  Inman  (1975),  Zohdy  (1975)  and  Pelton 
et  al.  (1978)  for  example.  In  the  RESIX  program,  an 
approximate  model  with  a  set  number  of  layers  is 
continually  modified  by  a  least-squares-error  pro¬ 
cedure,  providing  layer  thickness  and  resistivity 
values  that  best  match  the  field  data.  There  is  no 
assurance  that  the  model  is  unique,  only  that  it  is 
plausible  if  the  average  error  between  the  model 
curve  and  the  data  points  is  less  than  about  10%. 


RESULTS  AND  DISCUSSION 
New  England 

Kittery  and  Eastport,  Maine:  water  depth  studies 
Initial  observations  to  assess  equipment  per¬ 
formance  and  determine  the  influence  of  water 
depth  on  resistivity  values  were  made  at  Kittery 
and  Eastport,  Maine.  Profiling  was  done  in  New 
Haven  Harbor,  Connecticut,  since  our  results  could 
be  compared  to  the  large  amount  of  information 
available  on  depth  to  bedrock  (U.S.  Army  Corps  of 


Engineers  1981).  These  sites  were  also  selected 
because  the  seabed  electrical  properties  were 
similar  to  those  measured  in  some  subsea  perma¬ 
frost  zones  (Corwin  1983),  with  the  surface  of 
resistive  rock  simulating  the  top  of  ice-bonded 
permafrost. 

The  sites  at  Kittery  and  Eastport  were  chosen 
because  of  large  tidal  fluctuations  and  accessible 
outcrops  that  enabled  us  to  directly  measure  bed¬ 
rock  resistivity.  At  each  site  the  surface  sediments 
were  underlain  by  metamorphic  or  granitic  rocks 
that  outcropped  above  low  water.  No  information 
was  available  on  depth  to  bedrock  beneath  the 
cables,  but  the  many  outcrops  suggested  that 
bedrock  was  near  the  seabed.  Floating  and  bottom 
cables  were  secured  parallel  to  the  beach  and  data 
were  recorded  throughout  tidal  cycles  from  a  dock. 
Tidal  variations  were  about  7  m  at  Eastport  and 
about  2.5  m  at  Kittery  . 

The  resistivity  of  bedrock  outcrops  above  high 
water  and  the  resistivity  of  saturated  sediments 
and  bedrock  exposed  at  low  water  were  measured 
using  an  ac  magnetic  induction  technique  (the 
Geonics  EM-31)  (Arcone  et  al.  1979).  The  results 
are  summarized  in  Table  1 .  Tire  salinity  of  the  sea¬ 
water  was  about  30%o  and  the  resistivity  approx¬ 
imately  0.2  Q  m.  Preliminary  modeling  of  the 
Kittery  site  indicated  that,  for  an  array  with  a  spac- 
ings  up  to  60  m,  we  could  expect  apparent  resistiv¬ 
ity  to  vary  between  about  1  and  10  D  m  through¬ 
out  a  tidal  cycle. 

We  anticipated  that  water  depth  would  have  a 
significant  influence  on  observations  made  in  shal¬ 
low  coastal  waters.  Therefore,  the  response  of  a 
floating  array  at  the  Maine  sites  was  studied  both 
experimentally  and  by  using  a  mathematical  model 
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Apportnt  Rtiittmly  (£}  m) 


Table  1.  Resistivity  data  (Q  m)  from  shoreline  observations. 


Kitten i, 

Eastjvrt, 

Neti'  Haven, 

— . 

Maine 

Maine 

Connecticut 

Rock  outcrops  above  high  water 

— 

450  to  >10,000 

67  to  313 

(1-3  m) 

210  avg 

Rock  outcrops  exposed  at  low 

18  to  100 

45 

23  to  33 

water 

63  avg* 

28  avg 

Sediments  exposed  at  low 

2.4  to  4.1 

3.4  to  7.2 

— 

Water 

-.  ....  -  - 

_ 5.46  avg  .  _ 

_  ... 

*  Average  values  are  shown  where  more  than  10  measurements  were  made. 


Figure  3.  Model  illustrating  the  influ¬ 
ence  of  water  depth  on  observations  with 
a  floating  Wenner  array  for  the  Kitten/, 
Maine,  site.  The  resistivity  values  used 
were  0.25  Hi  nt  for  seawater,  2.0  Ti  m  for 
sand  and  20.0  Ti  m  for  the  rock.  The 
thickness  of  the  sand  layer  was  2.0  m. 


Figure  4.  Apparent  resistivity  data  obtained  with  a  floating 
Wenner  array  at  Kittery,  Maine  (solid  lines),  compared  with 
the  model  data  (dots)  shown  in  Figure  3. 


Figure  5.  Apparent  resistivity  data  obtained  with  a  floating  Wenner 
array  at  Eastport,  Maine. 
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Figure  6.  Influence  of  depth  of  seawater  on  apparent  resis¬ 
tivity.  Data  were  obtained  at  Eastport,  Maine,  with  a 
Wenner  array  on  the  seafloor. 


with  variables  being  water  depth  and  electrical 
properties  of  the  seafloor.  A  conservative  value  of 
20  Q  m  was  assigned  to  the  bedrock  so  as  to  exam¬ 
ine  our  limits  of  bedrock  sensitivity.  The  modeling 
results  are  shown  in  Figure  3,  with  data  from  field 
observations  made  during  the  tidal  cycles  shown 
in  Figures  4  and  5.  The  curves  all  have  the  same 
range  of  values  and  show  decreasing  sensitivity  to 
the  seabed  with  increasing  water  depth.  For  a  3-m 
water  depth,  an  a  spacing  of  at  least  25  m  is  re¬ 
quired  to  give  pa  values  above  the  resistivity  value 
of  the  sand.  Figure  4  shows  good  agreement  be¬ 
tween  field  data  and  model  at  the  Kittery  site  for 
low  and  high  water  limits. 

Similar  variations  in  resistivity  with  water  depth 
also  occur  in  data  collected  with  a  bottom  cable 
(Fig.  6).  In  addition  to  the  orderly  decrease  in 
sensitivity  to  the  seabed  with  increasing  water 
depth  seen  in  the  Kittery  data  and  model  and  in  the 
Eastport  data  for  the  floating  cable  (Fig.  3-5),  the 
bottom  cable  also  shows  a  noticeable  decrease  in 
sensitivity  to  water  depth  at  the  closely  spaced 
electrodes.  No  modeling  of  bottom  cable  responses 
was  attempted,  although  it  is  apparent  that  a  bot¬ 
tom  array  is  also  very  sensitive  to  water  depth 
(compare  a  =  10  m  in  Fig.  5  with  that  in  Fig.  6)  in 
areas  of  shallow  water.  Therefore,  on  the  basis  of 
these  results,  the  availability  of  models  for  inter¬ 
preting  data  only  from  surface  arrays,  and  the 
anticipated  damage  to  an  array  dragged  on  the 


seabed,  we  decided  to  use  a  floating  cable  for  our 
surveying  in  shallow  coastal  waters. 

New  Haven  Harbor:  bedrock  studies 

The  site  at  New  Haven  Harbor  was  selected 
because  of  the  substantial  information  available 
on  depth  of  bedrock  and  on  the  properties  and 
distribution  of  sediments.  The  information  is  based 
on  drilling  and  sampling  done  as  part  of  a  U.S. 
Army  Corps  of  Engineers  (1981)  study.  A  sample 
of  the  contoured  depth  to  bedrock  data  is  shown  in 
Figure  7. 

The  surveys  were  conducted  near  the  central 
part  of  the  harbor.  The  outcrop  resistivity  values 
(Table  1)  were  obtained  along  the  shore  and  are 
similar  to  the  values  for  the  Maine  sites.  The  elec¬ 
trical  properties  at  all  the  sites  are  similar  to  those 
expected  in  some  coastal  permafrost  settings,  par¬ 
ticularly  if  the  contact  between  ice-free  and  ice¬ 
bearing  materials  is  not  gradational. 

A  mathematical  model  (Mooney  et  al.  1966) 
was  first  used  to  investigate  the  theoretical  effect 
of  variations  in  depth  to  bedrock  for  resistivity  and 
depth  parameters  expected  in  New  Haven  Har¬ 
bor.  The  two  models  in  Figure  8  show  that  for  fixed 
water  depths  of  6  and  7  m,  a  variation  in  total  depth 
to  bedrock  of  47  m  will  result  in  only  a  1  Q  m 
change  for  a  surface  Wenner  array  with  a  =  50  m. 
This  fact,  coupled  with  the  large  influence  of  water 
depth  on  pa,  required  that  good  water  depth  con¬ 
trol  be  maintained  during  a  survey. 

Several  survey  lines  were  run  using  the  Wen¬ 
ner  marine  cable,  one  of  which  is  shown  in  Figure 
9.  This  profile  was  made  along  the  west  side  of  the 
channel  where  water  depth  was  fairly  constant  at 
6  ±  0.3  m.  Because  the  water  depth  variations  are 
small,  the  resistivity  profile  can  be  compared  di¬ 
rectly  with  observations  on  depth  to  rock  from  the 
Corps  of  Engineers  survey.  The  resulting  correla¬ 
tion  shown  in  Figure  9  is  very  good.  The  range  of 
pt  is  approximately  0.8  Q-m  for  a  total  depth  to 
bedrock  change  of  about  40  m,  which  agrees  with 
the  prediction  of  our  model  data  shown  in  Figure 
8. 

These  preliminary  observations  and  modeling 
demonstrated  the  utility  of  simple  dc  resistivity 
equipment  for  mapping  resistive  seabed  features 
in  shallow  coastal  waters.  The  low-power  Wenner 
system  provided  data  that  agreed  qualitatively 
with  variations  in  the  position  of  the  top  of  bed¬ 
rock  at  depths  as  great  as  50  m  below  sea  level.  The 
importance  of  having  good  control  on  water  depth 
in  shallow  water  surveys  was  also  illustrated. 

The  electrical  properties  of  these  sections  were 
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Prudhoe  Bay 

The  rates  of  coastal  erosion  and  associated 
modification  of  the  cold  terrestrial  permafrost  vary 
greatly  along  the  Beaufort  Sea  coastline  and  de¬ 
pend  on  material  type,  coastal  relief  and  exposure 
to  the  sea.  Three  sites  were  selected  for  this  study 
in  the  Prudhoe  Bay  area  and  are  shown  in  Figure 
10.  The  Prudhoe  Bay  area  has  low  coastal  erosion 
rates  compared  to  more  western  segments  of  the 
Beaufort  Sea,  where  rates  commonly  exceed  5  m/ 
yr  (Reimnitz  et  al.  1985).  Rates  in  Prudhoe  Bay  are 
closer  to  1-2  m/yr  and  are  probably  lower  than 
other  segments  of  the  coast  because  of  the  coarser 
material  and  lower  ground  ice  volumes  in  the 
upper  part  of  the  permafrost  section.  Infiltration  of 
sea  water  can  rapidly  influence  electrical  proper¬ 
ties  of  coarse-grained  material  if  voids  are  not  ice- 
filled.  In  general,  chemical  modification  by  salt 
should  be  much  less  rapid  in  fine-grained  than  in 
coarse-grained  sediment. 

The  soundings  on  land  were  all  made  asym¬ 
metrically  by  expanding  the  electrodes  in  one 
direction  parallel  to  the  water's  edge  at  a  spacings 
of  0.5, 1, 2, 3, 4, 5, 7, 10, 20, 30  and  50  m.  The  marine 
resistivity  data  were  gathered  every  25  m  along 
profiles  at  fixed  a  spacings  on  lines  normal  to  the 
coast.  The  offshore  a  spacings  commonly  used 
were  3,  6,  9,  12,  15,  30  and  46  m.  In  Figure  11 
contours  of  equal  values  of  apparent  resistivity  are 
mapped  as  a  function  of  distance  and  a  spacing. 
The  position  of  important  features  such  as  the 
coastal  bluff,  beach  and  shoreline  are  noted.  In¬ 
serts  with  expanded  scales  in  the  figures  illustrate 
the  rapid  changes  in  resistivity  seaward  of  the 
coastal  bluff.  No  insert  was  constructed  for  site  3 
since  changes  landward  were  gradual. 


Ditianct  (m) 


Figure  9.  Wenner  resistivity  profile  obtained  along  the  west  side  of  the  channel  in  New  Haven  Harbor  at  a  50- 
m  a  spacing. 


0  20  40  60  80 


Oepth  1o  Bedrock  (m) 

Figure  8.  Two  models  showing  the  range  in  resistivity 
values  expected  for  a  survey  in  New  Haven  Harbor  for  an 
a  spacing  of  50  m. 

within  the  range  of  those  expected  in  some  areas 
with  ice-bonded  subsea  permafrost.  In  some  cases 
ice-bonded  permafrost  will  be  more  variable  with 
gradational  contacts  and  layered  structure.  This 
naturally  will  make  detection  and  mapping  more 
difficult  for  the  dc  method,  as  it  would  for  any 
geophysical  technique.  These  positive  results  led 
us  to  the  permafrost  studies  reported  in  the  follow¬ 
ing  section. 


Figure  10.  Map  of  the  Prudhoe  Bay  area  showing  the  location  of  our  study  sites  and  the  orientation 
and  length  of  our  resistivity  profiles. 


Site  1 

This  site  (see  Fig.  1 0)  is  situated  several  hundred 
meters  east  of  the  West  Dock  in  Prudhoe  Bay.  It  has 
a  low  bluff  about  1.5  m  high  and  a  very  narrow  1- 
to  3-m  wide  beach.  The  fresh  nature  of  some  of  the 
exposures  suggests  active  erosion.  Only  a  small 
increase  in  sea  level  is  required  to  bring  waves  in 
direct  contact  with  the  bluff.  This  low  bluff  and  the 
adjacent  terrain  are  easily  overtopped  during 
storms,  as  indicated  by  the  large  amount  of  drift¬ 
wood  found  more  than  100  m  inland. 

Direct  observations  reported  by  Baker  (1987)  of 
the  depth  to  the  top  of  ice-bonded  permafrost 
offshore  along  one  control  line  are  shown  at  the 
top  of  Figure  11a.  This  line  provided  ideal  control 
for  the  marine  resistivity  study,  and  was  based  on 
drilling  and  probing.  Water  depth  data  included 
in  all  the  profiles  were  collected  during  our  resis¬ 
tivity  survey. 

Site  2 

This  site  is  about  3.7  km  northeast  of  the  East 
Dock  in  Prudhoe  Bay  (see  Fig.  10).  The  tallest 
coastal  bluffs  occur  here  (Fig.  lib,  top),  ranging 
from  2  to  3  m  in  height.  The  beach  is  about  6  m  wide 
and  is  covered  with  sod  blocks  that  slumped  from 
the  tundra  surface.  The  bluff  face  is  actively  erod¬ 


ing,  with  large  blocks  collapsing  over  fresh  expo¬ 
sures.  Wave  overtopping  of  the  bluff  appears  to  be 
very  limited,  with  only  minor  changes  in  the  vege¬ 
tation  patterns  near  the  bluff  crest  attributable  to 
salt  spray.  The  materials  exposed  in  the  bluff  are 
silty  sands,  which  seem  similar  to  the  exposed 
sediments  for  all  sites.  No  information  on  subsea 
permafrost  was  available  for  sites  2  and  3;  the 
subsurface  interpretations  shown  in  Figures  lib 
and  11c  are  the  results  of  computer  modeling  and 
will  be  discussed  below. 

Site  3 

This  site  is  approximately  2  km  east  of  the  West 
Dock  (see  Fig.  10).  This  location  is  different  from 
the  previous  sites  in  that  it  lacks  a  coastal  bluff  (Fig. 
11c,  top).  The  beach  formed  between  two  points  of 
land  and  grades  into  the  low-lying  adjacent  tun¬ 
dra.  Overtopping  of  the  beach  and  inland  flooding 
of  the  tundra  is  extensive  because  of  the  low  eleva¬ 
tion  of  the  tundra  surface.  Driftwood  can  be  found 
more  than  200  m  inland. 

Land  observations 

There  are  some  unique  features  in  the  on-land 
resistivity  data  that  appear  related  to  surface  re¬ 
lief,  beach  width  and  potential  for  surface  flood- 
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ing.  The  greatest  resistivities  were  observed  at  site  and  other  debris  indicates  is  the  inner  limit  of  most 

2  (Fig.  lib),  which  had  the  highest  bluffs  and  the  recent  flooding.  The  topography  in  this  area  sug- 

least  physical  indication  of  flooding  or  other  in-  gests  that  coastal  flooding  may  have  occurred  at 

land  disturbance.  A  maximum  apparent  resistiv-  this  site  for  many  years.  Values  greater  than  500 

ity  of  2800 11m  was  observed  65  m  inland  of  the  11m  are  not  found  at  depth  until  the  inner  limit  of 

bluff.  The  high  values  and  the  horizontal  disposi-  flooding  is  exceeded.  At  the  more  easterly  site  2, 

tion  of  the  contours  suggest  that  site  2  has  had  the  the  500-Hm  values  correspond  with  the  tall  bluffs 

least  modification.  In  contrast,  sites  1  and  3  show  and  higher-elevation  tundra  surface  that  mark  the 

the  influence  of  past  flooding  by  significantly  lower  inner  limit  of  flooding, 

values  of  near-surface  resistivities  and  the  sea¬ 
ward  sloping  of  the  resistivity  contours.  Near-  Seaward  obsewations 

surface  resistivities  at  site  3  are  less  than  6  Hm  at  There  are  some  common  patterns  in  all  data 

distances  more  than  175  m  inland  from  the  inner  from  seaward  of  the  coastal  bluffs.  No  apparent 

edge  of  the  active  beach,  a  distance  that  driftwood  resistivity  values  greater  than  100  Qm  were  ob- 
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a.  Site  1.  Contoured  data  for  an  expanded  segment  of  the  line  are  shown  in  the  lower 
part  of  the  figure.  The  position  of  the  top  of  ice-bonded  permafrost  (after  Baker  1 987) 
is  shown  in  the  upper  part  of  the  figure  as  a  solid  line.  The  computed  depth  to  the  top 
of  ice-bonded  permafrost  (shown  by  dots)  was  calculated  using  the  three-layer 
resistivity  models  at  the  top  of  the  figure  (resistivity  values  in  ohm-meters). 

Figure  1 1 .  Co)’.toured  apparent  resistivity  data  in  Prudhoc  Bay. 
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b.  Site  2.  Water  depth,  height  of  the  coastal  bluff  and  computed 
depth  to  the  top  of  the  ice-bonded  permafrost  are  shown  in  the 
cross  section  in  the  upper  part  of  the  figure.  Contoured  data  for 
an  expanded  segment  of  the  line  that  includes  the  beach  are 
shown  in  the  lower  part  of  the  figure. 

Figure  11  (cont'd).  Contoured  apparent  resistivity  data  in 
Prudhoe  Bay. 


served  at  even  the  largest  a  spacings.  The  values 
and  trends  were  similar  between  the  bluffs  and  the 
water  edge,  with  no  values  at  the  water  edge 
greater  than  50  Q  m  at  50-m  a  spacings;  most  val¬ 
ues  were  less  than  40  £lm. 

Variations  in  water  depth  along  the  study  lines 
were  usually  less  than  1  m,  the  total  depth  never 
exceeded  2  m,  and  the  water's  resistivity  was  al¬ 
ways  about  0.47  fl  m.  The  contoured  data  of  Fig- 
tire  11a  demonstrate  that  these  parameters  were 
not  sufficient  to  mask  resistivity  variations  at  depth. 


For  example,  the  larger  change  in  resistivity  at 
depth  in  Figure  11  occurs  in  a  zone  between  300 
and  400  m  where  there  is  less  than  a  0.5-m  vari¬ 
ation  in  water  depth.  This  small  change  in  water 
depth  cannot  account  for  the  more  than  6-fl  m 
resistivity  change  seen  across  the  zone  at  an  a 
spacing  of  46  m,  since  our  earlier  field  observa¬ 
tions  (Fig.  5)  and  theoretical  considerations  have 
shown  that  a  0.5-m  change  in  water  depth  would 
account  for  a  variation  of  less  than  1  £2  m.  The 
resistivity  changes  seen  at  the  small  a  spacings  (3-6 
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c.  Site  3.  Water  depth,  coastal  relief  and  computed  depth  to  the  top  of  ice-bonded 
permafrost  are  shown  in  the  upper  part  of  the  figure. 

Figure  11  (cant'd). 

m)  are  expected  for  the  small  changes  in  water  Apparent  resistivities  at  our  sites  fall  below  20 
depth  (Fig.  11a  and  c).  £2m  within  the  first  100  m  of  the  shore  at  46-m  a 

The  relatively  large  variations  in  resistivity  with  spacings. 

depth  and  the  configuration  of  the  contours  corre¬ 
spond  with  the  shape  of  the  top  of  ice-bonded  Modeling 

permafrost  observed  in  the  control  section  of  Fig-  Numerical  interpretation  of  the  offshore  ap- 

ure  11a.  However,  the  anomaly  in  the  resistivity  parent  resistivity  data  using  Baker's  (1987)  depth 
data  is  shoreward  of  the  break  (the  sharp  increase  profile  for  site  1  (Fig.  11a)  achieved  a  good  fit  with 

in  permafrost  depth  at  400  m)  in  the  control  data.  three-layer  models.  We  needed  3  three-layer 

This  discrepancy  may  be  caused  by  our  resistivity  models  to  interpret  the  entire  line;  they  are  shown, 

line  being  slightly  east  of  the  control  line,  but  is  along  with  the  calculated  position  of  the  top  of  ice- 

more  likely  attributable  to  limitations  in  resolu-  bonded  permafrost,  at  the  top  of  Figure  11a.  The 

tion  of  the  technique  itself,  which  are  discussed  first  layer  resistivity  was  held  constant  at  the 

below.  The  marine  resistivity  data  at  the  other  sites  measured  value  of  0.475  Qm.  Second  layer  resis- 

(Fig.  lib  and  c)  also  show  similar  breaks,  which  tivity  was  increased  with  water  depth  and  dis- 

are  taken  to  indicate  a  sudden  increase  in  depth  to  tance  from  shore.  A  resistivity  of  1 .1  Qm  was  used 

the  top  of  the  ice-bonded  permafrost.  The  offshore  when  water  depths  were  1 .5  m  and  less,  2.5  Q  m 

location  of  the  breaks  occurs  in  the  zone  where  for  the  1.6-  to  1.8-m  range,  and  2.7  Q  m  for  depths 

water  depths  are  approximately  1.5  m.  1.9  m  and  greater.  These  resistivities  seem  logical 

There  is  a  noticeable  gradient  in  all  the  resistiv-  because  of  the  possibility  of  salt  enrichment  of  the 

ity  data  that  fall  within  about  200  m  of  the  beach  bed  sediments  in  the  shallow  water  zone  when  salt 

and  appears  to  reflect  changes  in  the  electrical  is  rejected  during  formation  of  the  sea  ice.  At  great- 

properties  of  permafrost.  We  believe  this  to  be  er  water  depths  there  would  be  more  chance  for 

caused  by  warming  in  the  marine  environment.  mixing  and  freshening  of  the  bed  sediments. 
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The  only  significant  error  in  interpreted  depth 
is  at  the  400-m  station.  This  occurs  at  the  one  point 
where  the  permafrost  depth  drops  suddenly  and 
the  VVenner  array  would  have  spanned,  and  inte¬ 
grated  the  resistivity  over,  this  sudden  discontinu¬ 
ity  in  layering  depth.  The  greatest  change  in  depth 
occurs  over  a  50-m  distance,  which  approximately 
equals  the  largest  value  of  a  used  and  is  about  33% 
of  the  total  array  length.  At  all  other  stations, 
horizontal  changes  in  layering  are  sufficiently  small 
to  allow  this  one-dimensional  sounding  proce¬ 
dure  to  provide  good  resolution  for  this  two-di¬ 
mensional  cross  section. 

Apparent  resistivity  data  for  the  other  two  sites 
(2  and  3)  were  interpreted  using  the  model  resis¬ 
tivity  values  and  depth  relationships  developed 
for  site  1 .  The  calculated  position  of  the  top  of  ice- 
bonded  permafrost  is  shown  in  the  upper  parts  of 
Figures  11b  and  c.  These  models  are  not  greatly  in¬ 
fluenced  by  variations  in  resistivity  of  the  third 
layer.  Resistivity  of  the  ice-bonded  permafrost  can 
be  varied  from  200  to  1000  £2  m  without  significant 
change  in  the  calculated  depth  to  the  top  of  this 
layer.  However,  calculated  depths  of  the  ice- 
bonded  permafrost  are  extremely  sensitive  to 
changes  in  the  resistivity  of  the  second  layer.  For 
example,  where  the  observed  depth  to  the  top  of 
the  ice-bonded  permafrost  is  22  m  below  the  sea¬ 
bed,  a  ±0.5-£2  m  variation  in  the  resistivity  provid¬ 
ing  the  best  fit  can  cause  a  ±  6-m  variation  in  layer 
thickness.  Even  though  first-layer  parameters  are 
also  important,  they  are  not  a  problem  since  they 
can  be  directly  measured  at  the  time  of  a  survey. 

The  resistivity  values  obtained  for  the  frozen 
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Figure  12.  Model  resistivity  data  ( 50-m  a 
spacing)  for  various  water  depths  in  cases 
where  ice-bonded  permafrost  is  3-,  7-  and  20- 
m  depths  beneath  the  seabed. 


and  unfrozen  sediments  allow  us  to  evaluate  the 
ability  of  a  Wenner  array  with  a  maximum  elec¬ 
trode  spacing  of  50  m  to  perform  over  water  that 
exceeds  the  depths  encountered  here.  This  is  illus¬ 
trated  in  Figure  12,  which  shows  the  apparent  re¬ 
sistivity  calculated  for  a  50-m  a  spacing  as  a  func¬ 
tion  of  water  depth  in  cases  where  ice-bonded  per¬ 
mafrost  with  a  resistivity  of  900  £2  m  is  beneath 
three  sediment  (2.5  £2  m)  depths  of  3,  7  and  20  m. 
The  curves  show  that,  beyond  about  a  7-m  water 
depth,  there  is  little  sensitivity  to  the  total  depth  of 
permafrost.  This  was  not  a  problem  in  the  Prudhoe 
Bay  study  since  water  depths  did  not  exceed  3  m. 

SUMMARY  AND  CONCLUSIONS 

Using  dc  resistivity  data  from  an  arctic  coastal 
setting,  we  have  illustrated  the  impact  of  marine 
inundation  on  cold  permafrost  formed  on  land. 
The  floating  array  (a  =  50  m)  also  provided  a  tool 
for  measuring  the  electrical  properties  and  distri¬ 
bution  of  permafrost  in  shallow  coastal  waters 
where  ice-bonded  permafrost  may  not  be  more 
than  30  m  below  the  seabed  and  where  water 
depths  do  not  exceed  6  to  7  m.  The  technique  also 
has  application  to  subsea  bedrock  investigations 
under  certain  coastal  conditions.  Our  observa¬ 
tions  indicate  that  in  shallow  coastal  waters  (less 
than  about  7  m),  a  Wenner  seabed  cable  would 
offer  no  advantage  as  both  floating  and  seabed 
cables  are  about  equally  sensitive  to  changes  in 
water  depth.  However,  bottom  cables  should  have 
greater  sensitivity  to  seabed  changes  in  deeper 
water. 

Apparent  resistivity  observations  for  Wenner 
electrode  separations  up  to  50  m  made  along  sev¬ 
eral  cross  sections  extending  from  land  into  Prud¬ 
hoe  Bay  have  some  similar  characteristics:  1 )  resis¬ 
tivities  greater  than  20  £2  m  were  not  observed 
more  than  100  m  from  shore,  2)  maximum  resis¬ 
tivities  at  the  water's  edge  were  around  50  £2  m,  3) 
resistivities  greater  than  100  £2  m  were  not  found 
seaward  of  the  coastal  bluffs,  and  4)  resistivities 
just  inland  of  the  bluff  can  increase  above  the  sea¬ 
ward  values  by  an  order  of  magnitude.  Modeling 
indicates  that  real  resistivities  for  the  sediment 
over  ice-bonded  permafrost  ranged  from  1.1  to  2.7 
£2  m  and  that  the  permafrost  may  range  from  200 
to  1000  £2  m. 

Resistivity  contour  patterns  and  modeling  re¬ 
sults  correspond  well  with  the  drilling  and  pen¬ 
etrometer  observations  of  Baker  (1987);  thus,  this 
geophysical  approach  has  been  useful  in  this  geo- 


logical  setting  and  appears  to  have  applications 
for  scientific  and  engineering  investigations,  such 
as  offshore  pipeline  design  and  routing.  The  con¬ 
tours  and  modeling  suggest  that  a  rapid  increase 
in  the  depth  to  the  top  of  ice-bonded  permafrost 
seen  at  the  control  site  also  occurs  along  other 
study  lines  in  Prudhoe  Bay.  The  position  of  this 
zone  of  noticeable  increase  in  depth  to  the  top  of 
permafrost  seems  to  correspond  approximately 
with  the  1 .5-m  water  depth.  The  data  on  resistivity 
variations  within  200  to  300  m  of  the  coastline  may 
also  provide  useful  information  on  historical  ero¬ 
sion  rates. 
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